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We discuss aspects of the convex-roof extension of multipartite entanglement measures, that is, 
SL(2, <D) invariant tangles. We highlight two key concepts that contain valuable information about 
the tangle of a density matrix: the zero-polytope is a convex set of density matrices with vanishing 
tangle whereas the convex characteristic curve readily provides a non-trivial lower bound for the 
convex roof and serves as a tool for constructing the convex roof outside the zero-polytope. Both 
concepts are derived from the tangle for superpositions of the eigenstates of the density matrix. We 
illustrate their application by considering examples of density matrices for two-qubit and three-qubit 
states of rank 2, thereby pointing out both the power and the limitations of the concepts. 



I. INTRODUCTION 

One of the challenges in present quantum information 
theory is the quantification of multipartite entanglement 
both in pure and mixed states Once a measure for 
pure multipartite states is given (a "tangle" - that is, a 
polynomial invariant, see below), the corresponding mea- 
sure for mixed states can be obtained via the convex-roof 
extension Q. Although its definition is straightforward, 
the practical evaluation of the convex-roof extension is a 
difficult mathematical problem. To date a general ana- 
lytical method is known only for the concurrence of two- 
qubit mixed states [1, 0| ■ 

The convex roof is obtained by minimizing the average 
tangle of a given mixed state over all possible decompo- 
sitions of that state into pure states. As all elements in 
a decomposition are superpositions of the eigenstates of 
the given density matrix, the minimal tangle will depend 
on the possible tangle values for those superpositions. 

Recently, several authors have discussed how the en- 
tanglement of superpositions of some given pure states is 
related to the entanglement contained in those states [H, 
@, H B Hi- Here, we are considering a different prob- 
lem, namely how the entanglement of a mixed state de- 
pends on the entanglement of the superposition of its 
constituents. To this end, we introduce two concepts 
that facilitate the evaluation of mixed-state tangles for 
rank-2 states: the zero polytope and the convex charac- 
teristic curve. For a given mixed state p, the zero poly- 
tope is the convex set of those density matrices which 
have eigenstates in the range of p and, at the same time, 
have vanishing tangle. On the other hand, the convex 
characteristic curve is the largest convex function (of a 
suitable parameter) that does not exceed the tangle of 
any pure state in the range of p (for detailed explanation 
see below). 

The outline of this article is as follows. First, we intro- 
duce some terminology (Section II) and then give precise 
definitions of the zero polytope and the characteristic 
curve (Section III). There, we also discuss implications 
for the convex roof of the tangle in rank-2 states as well 



as the possible extension of these concepts to ranks larger 
than two. In Section IV we illustrate the application to 
examples of rank-2 states of two and three qubits. 

II. BASIC CONCEPTS 
A. Pure-state entanglement measures 

A function quantifying entanglement is required to be 
an entanglement monotone j 1 01 ] - Its crucial property is 
that it be non-increasing (on average) under stochastic 
local operations and classical communication (SLOCC). 
Stochastic local equivalence of two A-qubit states is rep- 
resented by convertibility under SL(2, C)® N transforma- 
tions [111, [l2j , and a homogeneous function on N qubits 
that is invariant under such operations is an entangle- 
ment monotone fl3l |. 

Here we investigate entanglement measures related 
to the concurrence introduced by Woottcrs and co- 
workers dQGl. 

The concurrence C(0ab) measures the degree of bipar- 
tite entanglement shared between the parties A and B in 
a pure two-qubit state \4*ab) G Ha ®Hb- In terms of 
the coefficients {0oo, 0oi> 0io> 4>n} of \4>ab) with respect 
to an orthonormal basis of product states it is defined as 

C(4>ab) = 2|0oo0n - 0oi0io | ■ (1) 

The concurrence is maximal for Bell states such as 
(| 00) + 1 1 1) )/ a/2 while it vanishes for factorized states. 

A measure for three-party entanglement has been in- 
troduced in Ref. [l5j]. The 3-tangle Tz{ip) of a three-qubit 
state ip £ Ha ®Hb <8> He can be expressed by using the 
wavefunction coefficients {"0oooj V'ooij • • • ; V'm} as 

r 3 = 4 \di -2d 2 + 4d 3 | 

d i = ^ooo^in + V'ooiV'iio + V'oioV'ioi + ^iooV&i 
0*2 = ^oooV'inV'oiiV'ioo + V'oooV'iiiV'ioiV'oio 

+^000^111^110^001 + V'OllV'100'0101^010 
+'0OllV'lOO'011O'0OOl + V'lOlV'OlO^llO^OOl 

di = V'oooV'iioV'ioi^oii + V'mV'ooiV'oioV'ioo ■ 
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For the GHZ state 

\GHZ) = ±= (|000) + |111)) (2) 

the 3-tangle becomes maximal: t 3 (GHZ) = 1. and it 
vanishes for any factorized state. It is most remark- 
able that there is a class of entangled three-qubit states 
for which T3 vanishes, i.e., these states are not SLOCC 
equivalent to any GHZ-type state [l2j ■ This class is rep- 
resented by the \ W) state 

\W) = -J= (|100) + |010) + |001)) , (3) 

and it is easy to check that t 3 (W) = 0. 

Both concurrence and the 3-tangle are the unique 
two- and three-qubit versions of more general polyno- 
mial measures invariant under local SL(2,C) transfor- 
mations [l3l . [l6j . These quantities are homogeneous of 
even degree in the coefficients of the wavefunction (with 
respect to a product basis) and play an important role as 
fundamental SL(2,C) invariants in invariant theory [13]. 
Therefore they serve as class-selective measures of (mul- 
tipartite) entanglement [l8|,[l||. In what follows, we refer 
to such polynomial invariants as tangles. 

B. Entanglement measures for mixed states — 
convex-roof extension 

Given a continuous real function g on the space of pure 
states n ex , this function can be extended to the space 
of mixed states in the following way [1, [2(3, [2l| . Let Q 
be the convex (and compact) set of normalized density 
operators. A state /) G Q can be written as a convex 
combination 

where the pure states it £ S! ox are extremal points of fi. 
The real function g u : Q — > K is a convex-roof extension 
of g : f2 ex — ► K if g u coincides with g on f2 ex and 

g u (p) := TDhx^pjgiiTj) . (5) 

where the minimum is taken over all decomposition of 
p. With this in mind it is formally straightforward to 
extend the definition of the concurrence to mixed two- 
qubit states p with a decomposition p = ^2Pjftf B into 
pure two-qubit states ir^ B of the parties A and B 

C{p) = min£>C(7r/ B ) . (6) 

Correspondingly, we have the 3-tangle for mixed three- 
qubit states p = J2Pj n f BC 

r 3 (p) = mm£>r 3 (7r/ BC ) (7) 

where itf BC denote pure three-qubit states. A decompo- 
sition for which the minimum of the respective function 
is realized, is called an optimal decomposition. 



III. ENTANGLEMENT OF SUPERPOSITIONS 
AND IMPLICATIONS FOR THE CONVEX ROOF 

It is in general difficult to evaluate the convex-roof ex- 
tension, i.e. the tangle of a given density matrix p. To 
this end, it is required to study the entanglement of pure 
states in the range of p, 1Z P , as can be seen from Defini- 
tion ([5]). Our starting point is the observation that the 
pure states \tp) p € 1Z P are superpositions of the eigen- 
states of p. 

A. Density matrices with vanishing tangle: the 
zero-polytope 

An important question is whether or not a tangle r 
vanishes for a given p. That is, we ask whether there 
is a decomposition of p into pure states {\4>j) } such 
that T(\tpj) ) = for all j = 0, . . . , r 2 — 1 where r is 
the rank of p. A polynomial SL(2, C) invariant r is a 
homogeneous function of the coefficients of \ip) p . The 
homogeneous degree of r (denoted by h) is an even pos- 
itive integer. Writing \ip} p = Y^j=o z j I'Aj)' * ms i m pli° s 
that t(|^i) ) is a polynomial of degree h in the complex 
variables Zo, ■ ■ ■ , z r -i. The situation is most easily un- 
derstood by considering rank-2 density matrices. Let us 
assume there is at least one eigenstate of p such that 
T (V'i) 7^ (otherwise we have trivially r(p) = 0). Since 
normalization does not play any role for this question, 
we then write 

\1)) p = \iM + z\^) . (8) 

The tangle is then a polynomial of degree h of a single 
complex variable z and consequently has precisely h zeros 
£1, . . . , (h- If both states have zero tangle we have to 
add the solution z = 00 to the h — 1 finite solutions of 
T(\ip) ) = 0. These zeros pif correspond to h pure states 

\Zi) p = |V'o)+Ci|^i) 

\z h ) p = IV'o) + Cfc IV'i) 

satisfying r(Zj) = for all j = 1, . . . , h. They span a 
polytope with h corners that contains precisely those den- 
sity matrices which can be decomposed into these states, 
i.e.. all convex combinations of ir 7 , . . . , ir 7 . All the den- 
sity matrices in this polytope have zero tangle (and out- 
side there is no other state with zero tangle) and therefore 
we call it the zero-polytope] its occurrence is generic for all 
polynomial SL(2, (D)-invariant entanglement monotones. 
This includes the concurrence, 3-tangle and multipartite 
entanglement monotones constructed in [l^] . 

We will now discuss briefly how the concept of the zero- 
polytope can be extended to higher-rank density matri- 
ces. As the tangle of homogeneous degree h evaluated 
on the pure state gives a polynomial of degree h in the 
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variables Zj, for z-x, ■ ■ ■ , z r —i kept fixed, we return to the 
situation explained above, and the remaining polynomial 
has precisely h zeros, Cj( z 2, ■ • ■ , Zr-i) with corresponding 
states 

r-i 

\Zj;Z2,...,Z r -l) = \lp )+Q{Z2,-.-,Z r -l)\lpl)+^2z j \'lp j ) ' 

k=1 

Upon varying z%, . . . , z r _i, each of these states will de- 
scribe an 2(r — 2) dimensional (real) manifold in the r 2 — 2 
dimensional manifold of states \tjj) (tp\. The convex hull 
of the union of these h manifolds contains exactly all den- 
sity matrices in 1Z P with zero tangle. 
Below we will discuss the zero-polytope for two explicit 
examples with rank r = 2. 

B. The convex characteristic curve 

It is well-known [22| that from the decomposition of a 
rank-r density matrix p = Y^j=o Pj l^j) (^j I its eigen- 
states {IV'j)}) an y other decomposition p = J2i \xi)(xi \ °f 
length m > r can be obtained with a unitary m x m ma- 
trix Uij via \xi) = YJj=o Uij y/pjH'j}- Therefore, the 
tangle of all pure states \ip) p 

W) P = X>i^> . En 2 = 1 

3=0 3 

can be regarded as a finger print of the tangle for all 
density matrices in 1Z P . This suggests to look for the 
minimum tangle of states \tp) p ; but rather than taking 
the global minimum, which is zero, we will look for local 
minima when the weights in the superposition are kept 
fixed. The degrees of freedom for this minimization are 
the relative phases in the superposition. In an appro- 
priate geometric representation, these minima describe 
characteristic curves of the tangle as a function of the 
weights. 

Let us discuss the simplest case first, that is, a rank- 2 
density matrix with its eigendecomposition 

p = Po\i>o)(i>o\+Pi\i>i)(H (9) 

with pj G IR, < pj < 1 such that po + Pi = 1. 
All pure states in 1Z P can be written as 

Mp = V«hto> + n/i^c^IVi) (10) 

where q, <p £ H and < <? < 1, < tp < 2ir. We may 
consider the tangle r of such a state (for a given phase <p) 
as a function of q in a r — q diagram (cf. also Fig. 1). This 
function represents a "characteristic curve" for r and is 
in general not convex. 

Clearly, each point (q, r(ip(q, ip))) in the r — q plane 
stands for the tangle of the pure state projector 
\ip(q, (p)) (ip(q, <p)\. Now consider an alternative length-2 
decomposition of p, p = p \(/>o) p (<Po\ +Pi\<f>i) p (<i>i\- 



The states \<fro) and \<j>\) p have the we weights q , 
qi and the relative phases tpo and (fx- As these states 
form a decomposition of p their parameters obey po = 
Polo + PiQi- Moreover, the phases need to be adjusted 
properly. Note that, for arbitrary phases, a density ma- 
trix po |V>o} p (^o| +Pi H'i) p (ipi\ +a\ipo) p (ipi\ +h.c. is ob- 
tained which, in the eigenbasis, has the correct diagonal 
elements. 

The average tangle of the decomposition p = 
Po \<l>o) p ((t>o \ + Pi \<f>i) p {4>i\ is visualized in the r - q plane 
by connecting the points (q , T{q ,tp )) and (qi, r(qi, tpi)) 
corresponding to the tangle of the states \<fio) p , \<t>i} p by 
a straight line. 

Its value po r (0o) +Pit{4>i) is assumed for the abscissa 
value q = po = po9o +P1Q1 ■ The generalization to decom- 
positions of p of length > 2 is straightforward. Hence, a 
given decomposition of p can be assigned a point in the 
t — q plane with abscissa po and ordinate equal to the 
average tangle r that results from convexly combining 
the tangles of the 

pure states in the decomposition. 




FIG. 1: Characteristic curves r(q, ip) in the r — q plane (solid 
gray lines) for the superpositions in Eq. (llOp . The (blue) 
dashed line represents the minimum f(q) from Eq. and 
the solid (red) line its function convex hull r*(q), i.e., the con- 
vex characteristic curve (for details see text). 
The point with label p indicates the average tangle of a par- 
ticular length-2 decomposition of the density matrix p — 
Po \<f>o) + Pi 1 0i ) - Moreover we illustrate the actual tangle of 
this density matrix r(p) and its lower bound r*(po) according 
to the convex characteristic curve. 

The actual value of the tangle r(p) is the minimum for 
all possible decompositions. In our geometric visualiza- 
tion this means to find the smallest possible r(po) that 
may result from a convex combination of at most r 2 pure 
states. In order to find this value of r (or at least ap- 
proximate it) one may try the following approach: first 
minimize the pure-state tangle r(q, tp) for each argument 

q 

f(q) = mm[r(il>(q,<p))] (11) 

v 

and then try to find decompositions from these minimal 
values f(q). That is, we need to investigate convex com- 
binations of pure states and their corresponding tangle. 
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However, we note that the function f(q) is in general not 
convex. 

Therefore we define r*(q) as the function convex hull 
of f{q) (i.e., the largest convex function that does not ex- 
ceed r{q,tp)) and call it the convex characteristic curve. 
It is not difficult to see that the convex characteristic 
curve provides a non-trivial lower bound for the tangle 
r{p) of the initial rank-2 density matrix p. Any (nor- 
malized) vector |</>° pt ) p in an optimal decomposition of 
p = ^2p° pt | < / > j Pt ) p (0j Pt | is characterized by some qj and 
a phase tpj in analogy with Eq. (|10p. A lower bound 
for the tangle of the pure state |</>j Pt ) (0j Pt | i s given by 

the convex characteristic curve: r(0° pt ) > T*(qj). Conse- 
quently, the tangle of p is bounded from below by the con- 
vex combination Pj Pt T * (ij) °f ^ ne points {(Qj,T*(qj))} 
on the convex characteristic curve. As the weight of l^o) 
in p equals po we find 

<p) = E^r^D ^ E^ PV fe) > t*(po) ■ (12) 



This is a tight estimate of r(p). The equal sign in the 
first inequality in Eq. (TT2]) holds if and only if there is an 
optimal decomposition such that all points correspond- 
ing to the states of this decomposition lie on the convex 
characteristic curve. For the equal sign in the second 
inequality in Eq. (|12[) . cither r* must be affine or all q{ 
have to coincide. 

If the rank of the given density matrix is larger than 
two, the convex characteristic curve can be generalized to 
a convex manifold. Suppose we have a family of rank-r 
density matrices with its eigendecomposition 



P = Po|V'o}(V>o| + ---+Pr-l|VV-l}(VV-l| (13) 



and Y^=o Pi = 1 ■ Then we define 



f(q ,...,q r - 2 ) = mmT[ilj(qo,...,q r -2;<Pi,---,<Pr-l) 

r-1 

\ipiq ,...,q r - 2 ;ipi,...,(p r -i)) = Vqo H'o) + E V<h ellf3 ' 1^') 

3=1 



(14) 
(15) 



r 



with qj G IR such that X^=o Qj = !■ The convex char- but not necessarily a tight one 
acteristic manifold r*(q , . . . , q r - 2 ) is then defined as the 
function convex hull of f(qo, . . . , q r — 2). 

In analogy with the rank-2 case, the value of 
t*(pq, . . . ,p r _2) gives a strict lower bound for the tan- 
gle of the density matrix in Eq. (fT3|) . However, due to 
its higher dimensionality it is much more difficult than in 
the rank-2 case to completely characterize this manifold 
and use it for estimates. Even as a tool for a numerical 
or graphical method it will be rather hard to handle. 

In order to illustrate the relevance of the convex char- 
acteristic curves /manifolds, but to also highlight caveats, 
we discuss two examples on specific rank-2 density ma- 
trices for two and three qubits. 



A. Mixed three-qubit states with GHZ and W 
components 

For three qubits, the GHZ state and the W state are 
given by Eqs. ([2]) and ([3]), respectively. We note that 
the 3-tanglc of the GHZ state is maximal: t^(GHZ) = 1 
while it vanishes for the W state t 3 (W) = More- 
over, we have (GHZ\W) = 0. 

In this example, we consider the family of mixed states 



p(p) = p 7T GH Z + (l-p)7T W , 



(16) 



IV. EXAMPLES 

In the following we illustrate the concepts introduced 
above in two examples: three-qubit mixtures of GHZ and 
W states [23|, and a generic two-qubit example. While in 
the highly symmetric three-qubit case the convex char- 
acteristic curve leads directly to an analytic solution for 
the convex roof of a whole family of states, the two-qubit 
example shows that in general the convex characteristic 
curve will give a lower bound for the mixed-state tangle 



where p is a real number and < p < 1. Our goal is to 
find the 3-tangle r^ipip)) for all values of p. 

According to our discussion above, all elements of the 
optimal decompositions of pip) are linear combinations 
of \GHZ) and \W) 



\Z(j>,<p)) = Vp\GHZ) - e^y/T^p\W) . (17) 

Let us first find the corners of the zero-polytope. By 
using the formula for the 3-tangle Eq. ([2]) for the states 
\Z(p, ip)) we obtain 



r 3 (Z(p, (fi)) 



8V$ 



pfe Mv 



. (18) 
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This equation is obeyed by the W state (p = 0) 
as by the three states 



Z 3 n )=y/to\GHZ) 



- e^^l - ~p~ \W) 



with j = 0, 1, 2 and 



Po 



4^2 
3 + 4\/2 



= 0.626851. 



well 



(19) 



(20) 



These four states define the zero-polytopc Sq for the fam- 
ily p(p) (which in this case is a simplex; see also Fig. [3]). 
We note that all states inside Sq have optimal decomposi- 
tions of length four while the states on the surface of So 
have three- vector optimal decompositions. It is worth- 
while mentioning that we have found many more mixed 
states with vanishing 3-tangle than just those that be- 
long to the family p(p). Indeed, as explained above, Sq 
contains all density matrices in 1Z P with zero 3-tangle. 

Interestingly, it is possible to determine the complete 
convex roof for p{p) as we will show now. To this end, 
we need to consider the corresponding convex character- 
istic curve. We have plotted T${Z{p, p)) in Fig. [2] and its 




FIG. 2: The 3-tangle of the states \Z(p,<p)) in Eq. JT7J for 
various values of <p = 7 • ^ (from top to bottom: 7 = 
1/2,1/3,1/5,1/10,0). 



minimum (solid black line for tp = 0) 
f 3 (p) = mmT 3 {Z(p,p)) 

, 8^ 



Vp( 1 -p) 3 



(21) 
(22) 



The function f 3 (p) is not convex for < p < po and 
0.825... < p < 1. Its function convex hull, i.e., the convex 
characteristic curve for p(p) in Eq. (fl"6|) is given by 



r 3 * (P) 



where 



, 
hip) 



0<p< Po 
Po<P<Pi 



Pi 



l-(l-p)(| + iV465) , pi < P < 1 

(23) 

~1 = 0.70868... (24) 



As we have explained in the preceding section r*(p) 
represents a lower bound for r{p{p)), and it coincides 
with the 3-tangle for those values p for which the corre- 
sponding T*(p) is realized by at least one decomposition. 
For the states p(p) it is not difficult to specify a decom- 
position with an average 3-tangle equal to t* (p) (see also 
Ref. [H). 

In the interval po < p < p\ we have the "star-shaped" 
three- vector decomposition Z 



Pip) 



1 



(\Z a )(Z°\ + \Z 1 )(Z 1 \ + \Z 2 )(Z 2 \)(25) 



\z°) = Vp\ghz) - ^T^p~\w) 
\z 2 ) 



^\GHZ) - e^v/l - p\W) 
^\GHZ) - e^^l~^p\W) 



(26) 



For p < po a four- vector decomposition is optimal that 
combines the decomposition Z in Eq. (f2"5|) for po with 
the W state: 



Pip) = — PiPo) 
Po 



1-P- 

Po 



7T\V 



(27) 



Finally, for po < p < 1, we need to combine the decom- 
position Z for pi with the GHZ state: 



Pip) 



I-P 
1-Pi 



P(Pi) 



P-Pi 
I -Pi 



ttqhz 



(28) 



We stress that, as the convex characteristic curve T 3 {p) 
gives a strict lower bound for T 3 {p{p))^ the existence of de- 
compositions attaining these values means that we have 
found the convex roof analytically. It is worth noticing 
that the above given optimal decomposition is unique. 



GHZ (p=l) 




Pi 
Po 



310 



FIG. 3: Bloch sphere for the two-dimensional space spanned 
by the GHZ state and the W state. The simplex So contains 
all states with vanishing 3-tangle. The "leaves" between po 
and pi represent sets of constant 3-tangle. In the simplex Si , 
the 3-tangle is an afhne function (for further explanation, see 
text). 

The extension of the analytic convex roof to the en- 
tire simplices spanned by the decomposition vectors (see 
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Fig [3]) follows from the theorem that the vectors in an 
optimal decomposition are also optimal for the entire 
simplex they span (this applies also to states forming 
a subset of an optimal decomposition) 0). 

B. Generic example for a rank-2 mixed two-qubit 
state 

In the previous example we have seen that studying 
the convex characteristic curve may lead to a complete 
solution for the convex roof of a family of mixed states. In 
general, however, this may not be expected. The reason 
why the characteristic curve will typically fail to give the 
exact convex roof is because the phases ipj for which the 
minimum in Eqs. (fTTjl . (fT4"| is reached, do not necessarily 
admit a decomposition of p. 

Nevertheless one obtains valuable information about a 
given mixed state. In order to illustrate this, we describe 
a rather generic two-qubit case. Let us consider the states 

\I) = -L (2 100) + |11)) (29) 

= -L (|00) + |01) -2 |11)) (30) 
and their convex combinations 

02 (?) = V 7TJ + (1-P)7T// , (31) 

We note that (I\II) = 0. This example can easily be 
treated analytically by using Wootters' method [|[. 

In particular, we find vanishing concurrence for po = 
5/11 (sec Fig. gl). 




FIG. 4: Solid (blue) line: exact solution for the concur- 
rence C(p2(p)) according to Eqs. (|29l) - (|31[) . For comparison: 
the convex characteristic curve of pi{p) (dashed red line; cf. 
also Fig. 5). The exact concurrence C(p(p)) has a zero for 
Po = 5/11. 

Now let us first evaluate the zero-polytope (which is a 
line in the Bloch sphere in this case). The equation 

C( y /p\I) + e i *y/T=]>\II))=0 (32) 

has two real solutions for p: pi — 5/29 with (p — 7r and 
Pi = 10/13 with tp = 0. ^From this we can conclude im- 
mediately that there is exactly one value of p for which 



the mixed-state concurrence of p2(p) vanishes: it is the 
intersection point of the zero-polytope with the line in the 
Bloch sphere that represents the family pi{p)- Straight- 
forward algebra shows that this intersection occurs at 
Po = 5/11. 

We see that the zero-polytope indeed gives us exact 
information about a subset of the family piip) only (here: 
about a single element). Note however that this subset 
can even be empty in the most general case, when e.g. 
the wave function coefficients of ([2"9"| have a non-zero 
relative phase modulo n. 




0.2 0.4 0.6 0.8 1 
P 



FIG. 5: Concurrence for superpositions of states |7) and 
1 77) (solid gray lines) and convex characteristic curve C*(p) 
(dashed red line). The convex characteristic curve coincides 
with the abscissa (i.e., C*(p) = 0) between the two zeros of 
Eq. (02J, pi = 5/29 and p 2 = 10/13. 

In Fig. 5 we show how the convex characteristic curve 
of p2 (p) arises from the concurrence of superpositions of 
the states \I) and \II). Analytically, the convex charac- 
teristic curve reads: 

r c-(p) , q<p<pi 

C*( P ) = { , p 1 <p<p 2 (33) 

{ C+(p) , ft <p< 1 

where C±(p) = C(y/p\I) ±y/T=p\II)). Clearly, C*(p) 
must vanish for all values of p between the zeros accord- 
ing to Eq. ([32]). We observe (cf. Fig. 4) that C*{p) is 
indeed a lower bound to the exact concurrence. How- 
ever, the latter is strongly underestimated by the convex 
characteristic curve for almost all p values. 

V. DISCUSSION 

We have introduced two concepts that serve to facil- 
itate the investigation of the convex roof of tangles for 
multipartite mixed qubit states. Both concepts are based 
on the tangle of pure-state superpositions, that is, of 
those pure states that lie in the range of the original 
density matrix. 

The zero-polytope is formed by the convex set of den- 
sity matrices with vanishing tangle in the range of the 
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original mixed state under consideration. It can always 
be determined exactly. Then, the remaining task is the 
calculation of the convex roof for all remaining density 
matrices with non-zero tangle. 

To this end, we have proposed the concept of the con- 
vex characteristic curve (or manifold) . It leads to a non- 
trivial and possibly tight lower bound for the tangle of 
a given density matrix. In particular, it coincides with 
the convex roof if there is a decomposition of the original 
density matrix into pure states whose tangles lie on the 
convex characteristic curve. 

In order to illustrate the power of these concepts we 
have discussed two simple examples - three-qubit mixed 
states with GHZ and W components [23| , and a generic 
rank-2 two-qubit state. For the GHZ/W mixtures the 
convex characteristic curve shows analytically i) that the 
results of Ref. [23[ indeed represent the convex roof for 
that family of states, and ii) provide an example where 
the convex characteristic curve provides the complete so- 
lution for a convex roof - namely, due to the fact that 
for each point of the convex characteristic curve there ex- 
ists a decomposition that attains this lower bound. On 
the other hand, the two-qubit example displays that in a 



more general case, one may obtain non-trivial estimates 
for the convex roof. However, as in that case there hardly 
exists an optimal decomposition with all its elements on 
the convex characteristic curve, on must not expect that 
the obtained estimates for the convex roof are close to 
the exact value. 

While, in principle, both concepts can be generalized 
to mixed states of arbitrary rank, their application is 
particularly adapted to the case of rank-2 states. 
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